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Abstract 

This paper is devoted to investigate the exact solutions of Bianchi 
type I spacetime in the context of f{R,T) gravity [T], where f{R,T) 
is an arbitrary function of Ricci scalar R and trace of the energy 
momentum tensor T. For this purpose, we find two exact solutions 
using the assumption of constant deceleration parameter and variation 
law of Hubble parameter. The obtained solutions correspond to two 
different models of the universe. The physical behavior of these models 
is also discussed. 
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1 Introduction 

The most popular issue in the modern day cosmology is the current expansion 
of universe. It is now evident from observational and theoretical facts that our 
universe is in the phase of accelerated expansion Bl-lII]- The phenomenon 
of dark energy and dark matter is another topic of discussion |12]-[l9]. It was 
Einstein who first gave the concept of dark energy and introduced the small 
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positive cosmological constant. But after sometime, he remarked it as the 
biggest mistake in his life. However, it is now thought that the cosmological 
constant may be a suitable candidate for dark energy. Another proposal to 
justify the current expansion of universe comes from modihed or alternative 
theories of gravity. /(T) theory of gravity is one such example which has 
been recently developed. This theory is a generalized version of teleparal¬ 
lel gravity in which Weitzenbock connection is used instead of Levi-Civita 
connection. The interesting feature of the theory is that it may explain the 
current acceleration without involving dark energy. A considerable amount 
of work has been done in this theory so far [20] . Another interesting modified 
theory is f{R) theory of gravity which involves a general function of Ricci 
scalar in standard Einstein-Hilbert Lagrangian. Some review articles [21] can 
be helpful to understand the theory. 

Many authors have investigated f{R) gravity in different contexts |22j - 
[35] . Spherically symmetric solutions are most commonly studied solutions 
due to their closeness to the nature. Multamaki and Vilja [36] explored vac¬ 
uum and perfect fluid solutions of spherically symmetric spacetime in metric 
version of this theory. They used the assumption of constant scalar curvature 
and found that the solutions corresponded to the already existing solutions in 
general relativity (GR). Noether symmetries have been used by Capozziello 
et ah [37] to study spherically symmetric solutions in f{R) gravity. Simi¬ 
larly many interesting results have been found using spherical symmetry in 
f{R) gravity [3B]. Cylindrically symmetric vacuum and non-vacuum solu¬ 
tions have also been explored in this theory [39]. Sharif and Shamir [ID] 
found plane symmetric solutions. The same authors [H] discussed the so¬ 
lutions of Bianchi types I and V cosmologies for vacuum and non-vacuum 
cases. Conserved quantities in f{R) gravity via Noether symmetry approach 
have been recently calculated [12]. 

In a recent paper [1], Harko et al. proposed a new generalized theory 
known as f{R,T) gravity. In this theory, gravitational Lagrangian involves 
an arbitrary function of the scalar curvature R and the trace of the en¬ 
ergy momentum tensor T. Myrzakulov [H] discussed f{R,T) gravity in 
which he explicitly presented point like Lagrangians. Sharif and Zubair [H] 
studied the laws of thermodynamics in this theory. The same authors [15] 
investigated holographic and agegraphic f{R,T) models. Houndjo [ID] re¬ 
constructed f{R,T) gravity by taking f{R,T) = fi{R) + f 2 {T) and it was 
proved that f{R,T) gravity allowed transition of matter from dominated 
phase to an acceleration phase. Thus it is hoped that f{R,T) gravity may 
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explain the resent phase of cosmic acceleration of our universe. This the¬ 
ory can be used to explore many issues and may provide some satisfactory 
results. 

The isotropic models are considered to be most suitable to study large 
scale structure of the universe. However, it is believed that the early uni¬ 
verse may not have been exactly uniform. This prediction motivates us to 
describe the early stages of the universe with the models having anisotropic 
background. Thus, the existence of anisotropy in early phases of the universe 
is an interesting phenomenon to investigate. A Bianchi Type I cosmological 
model, being the generalization of flat Friedmann-Robertson-Walker (FRW) 
model, is one of the simplest models of the anisotropic universe. Therefore, 
it seems interesting to explore Bianchi type models in the context of /(R, T) 
gravity. Adhav m investigated the exact solutions of f{R,T) held equa¬ 
tions for locally rotationally symmetric Bianchi type I spacetime. Reddy 
et ah IH] explored the solutions of Bianchi type III spacetime using the 
law of variation of Hubble’s parameter. Bianchi type III dark energy model 
in the presence of perfect huid source has been reported [IH]. Ahmed and 
Pradhan [SO] studied Bianchi Type V cosmology in this theory by involving 
the cosmological constant in the held equations. Naidu et ah [5T] gave the 
solutions of Bianchi type V bulk viscous string cosmological model. 

In this paper, we are focussed to investigate the exact solutions of Bianchi 
type I spacetime in the framework of f{R,T) gravity. The plan of paper is 
as follows: In section 2 , we give some basics of f{R,T) gravity. Section 
3 provides the exact solutions for Bianchi type I spacetime. Concluding 
remarks are given in the last section. 


2 Some Basics of /(i?, T) Gravity 


The f{R,T) theory of gravity is the generalization or modihcation of GR. 
The action for this theory is given by |1] 



( 1 ) 


where f{R, T) is an arbitrary function of the Ricci scalar R and the trace T 
of energy momentum tensor while Lm is the usual matter Lagrangian. It 
is worth mentioning that if we replace f{R,T) with f{R), we get the action 
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for f{R) gravity and replacement of f{R,T) with R leads to the action of 
GR. The energy momentnm tensor is defined as [S2] 


Here we assnme that the dependance of matter Lagrangian is merely on the 
metric tensor g^i, rather than its derivatives. In this case, we get 

T^u = LmQfiu - (3) 


The f{R,T) gravity held eqnations are obtained by varying the action S in 
Eq.(IT]) with respect to the metric tensor 


fn{R, T)R,,~f{R, T)g,,-{V,V,-g,.n)fn{R, T) 

where denotes the covariant derivative and 


/t^R) T){Tfj_u+Q, 

(4) 


□ = V^V^, /k(R,T) 


dfR{R,T) 

dR 


fAR,T) 


dfR{R,T) 


Contraction of (|1]) yields 


/«(R, T)R + 3n/B(R, T) - 2f{R, T) = rT - MR, T){T + 0), (5) 


where 0 = This is an important eqnation becanse it provides a rela¬ 

tionship between Ricci scalar R and the trace T of energy momentnm tensor. 
Using matter Lagrangian L^, the standard matter energy-momentnm tensor 
is derived as 

M = (p + p)u^Uu - pg^u, (6) 

where 0, 0, 0) is the fonr-velocity in co-moving coordinates and 

p and p denote energy density and pressure of the huid respectively. Perfect 
huids problems involving energy density and pressure are not any easy task 
to deal with. Moreover, there does not exist any unique dehnition for matter 
Lagrangian. Thus we can assume the matter Lagrangian as = —p which 
gives 

= -pPfMu - 2T^i., (7) 
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and consequently the field equations (jl]) take the form 


( 8 ) 

It is mentioned here that these field equations depend on the physical nature 
of matter field. Many theoretical models corresponding to different matter 
contributions for f{R, T) gravity are possible. However, Harko et ah [1] gave 
three classes of these models 

r R+2fiT), 

f{R,T)={ MR) + MT), 

I U{R) + UR)UT). 


In this paper we are focussed to the first class, i.e. f{R,T) = R + 2/(T). 
For this model the field equations become 


Rfiu — -^Rdiiu — + 2/'(T)T^j, + 


f{T) + 2pf{T) 




(9) 


where prime represents derivative with respect to T. 


3 Exact Solutions of Bianchi Type / Universe 

In this section, we shall find exact solutions of Bianchi I spacetime in /(i?, T) 
gravity. For the sake of simplicity, we use natural system of units {G = c = 
1) and /(T) = AT, where A is an arbitrary constant. For Bianchi type I 
spacetime, the line element is given by 

ds^ = dt^ — A^{t)dx^ — B‘^{t)dy^ — C^{t)dz‘^, ( 10 ) 


where A, B and C are defined as cosmic scale factors. The Bianchi I Ricci 
scalar turns out to be 


R 


-2 


A ^ B ^ C ^ Ab ^ BC ^ CA 


( 11 ) 


where dot denotes derivative with respect to t. 
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Using Eq.(|n]), we get four independent field equations, 


AB BC ^ CA 
'AB^Wj^'CA 
B C BC 

C A CA 
c ^ I ^ lu 
A B Ab 
A^B^AB 


(Stt + 3A)p — Xp, 
Xp — (Stt + 3A)p, 
Xp — (Stt + 3A)p, 
Xp — (Stt + 3A)p. 


( 12 ) 

(13) 

(14) 
(16) 


These are four non-linear differential equations with Eve unknowns namely 
A, B, C, p and p. Subtracting Eq. ffTTl) from Eq. ffT^ . Eq. ffT^ from Eq. ffTT]) 
and Eq. ffT^ from Eq. ffT^ . we get respectively 


A 

A 

B 

B 

A 

A 


B C/A 
b^c[a 
C AfB 
C^A\B 
C 5/i 
C^B\A 



(16) 

(17) 

(18) 


These equations imply that 


— — di Gxp 

A 

— = d 2 exp 
±) 

— = ds exp 


f dt 
f dt 

(U < 

f dt 

C3 / ^ 


(19) 

( 20 ) 
( 21 ) 


where Ci, C 2 , C 3 and di, d 2 , d^ are integration constants which satisfy the 
following relation 

Cl + C 2 + C 3 = 0 , did2d3 = 1 . ( 22 ) 

Using Eqs. (ll9h - (|2T1) . we can write the unknown metric functions in an explicit 
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way 


A = api exp gi 


B = ap 2 exp q 2 


C = aps exp qs 


where 


Pi = {di d2 )3, p2 = {did2 )3, p^ = {did2)^ 


2ci + C 2 Cl — C 2 Cl + 2 c2 , 

Qi = -^, 92 = — ^ —, 93 =-^-• (27 

It is mentioned here that pi, p 2 , Ps and qi, q 2 , qs also satisfy the relation 

P 1 P 2 P 3 = 1, 9i + 92 + 93 = 0. (28 


3.1 Some Important Physical Parameters 

Now we present some important dehnitions of physical parameters. The 
average scale factor a and volume scale factor V are dehned as 

a = ^ABC, V = a^ = ABC. (29) 

The generalized mean Hubble parameter H is given by 

H =^-{H^ + H2 + H^), (30) 

where TTi = ^, H 2 = ^ are dehned as the directional Hubble 

parameters in the directions of x, y and axis respectively. The mean 
anisotropy parameter A is 

i=l ^ ' 

The expansion scalar Q and shear scalar are dehned as follows 


2 3 \a) 


^ V c / AB BC CA ^ 
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where 


(34) 


with hfj,y = Qfj^y — u^Uy defiiied as the projection tensor. The deceleration pa¬ 
rameter q is the measnre of the cosmic accelerated expansion of the nniverse. 
It is defined as 

da 

^ = ~~ 

The behavior of the nniverse models is determined by the sign of q. The 
positive valne of deceleration parameter snggests a decelerating model while 
the negative valne indicates inflation. Since there are fonr eqnations flT^ - 
03 and five nnknowns, so we need an additional constraint to solve them. 
Here we use a well-known relation [53] between the average scale factor a 
and Hubble parameter H to solve the equations, 



H = la-^, 


(36) 


where I and n are positive constants. This is an important relation because it 
yields a constant value of deceleration parameter and consequently we obtain 
power law and exponential models of universe. Using Eqs. fl30|l and fl36|l . we 
get 

d = la^-^ (37) 

and the deceleration parameter becomes 


q = n — 1. 

Integrating Eq. (l57D . it follows that 

a = (nit + ki)^, n ^ 0, 

and 

a = k 2 exp(lt), n = 0, 
where ki and ^2 are constants of integration. 


(38) 


(39) 

(40) 


3.2 Singular Model of the Universe 

Here we investigate the model of universe when n 7 ^ 0, i.e., a = {nit + /ci)’ 
In this case, the metric coefficients A, B and C take the form 


A = pi{nlt + ki)^ exp 
B = p 2 {nlt + ki)^ exp 
C = p^{nlt + ki)^ exp 


qi{nlt + ki)' 
l{n — 3) 

q2{nlt + ki)' 
l{n — 3) 

q^{nlt + ki)' 


77- 7 ^ 3, 
n 7 ^ 3, 
n 7 ^ 3. 


l{n — 3) 

The directional Hubble parameters ifj (i = 1, 2, 3) turn out to be 

H = ^ 

nit + ki {nit + ki)n 

The mean generalized Hubble parameter and volume scale factor are 

I 


H = 


V = {nit + fci)". 


nit + ki 

The mean anisotropy parameter becomes 

+ qs" 


A = 


3 P{nlt + ki)(^ 2n)/n 
The expansion scalar and shear scalar for this model are given by 


9 = 


31 o + q2^ + 


a2 = 


(41) 

(42) 

(43) 


(44) 


(45) 


(46) 


nlt + ki “ 2 {nlt + kiY/'^ ' 

Using Eqs. the energy density of the universe turns out to be 


(47) 


P = 


12(A + 27r)(A + 47r) _ 
- A 


4(A + dvr) 


3/2 


{nit + fci)2 


+ 


gig2 + g2g3 + q^qi 


{nit + /ci)' 


3/2(1-n) ^q,^ + q,^ + q,^ 


{nit + ki 


{nit + ki 


(48) 
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Figure 1: Behavior of energy density and pressure versus time for t > 0 with 
n = 2, A = 1, / = 1, ki = 0, qi = 1 = q 2 and q^ = — 2 . 


while the pressure of the universe becomes 


P = 


12(A + 27r)(A + 47r) 

3/2(1 -n) 


Att 


(3A + Stt) 


{nit + /ci )2 


+ 


3/2 ^ qiq2 + q2(l3 + QsQi 

{nit + ki)"^ {nit + ki)^ 
qi^ + + gsM ' 

{nit + ki)n J 


(49) 


The plots of p, p and equation of state parameter w = p/p against time 
coordinate t are shown in hgure 1 and 2 respectively. It is evident from 
hgure 2 that w —| as / —)■ cxo. Thus the model corresponds to a radiation 
dominated universe as the time grows. 


3.3 Non-singular Model of the Universe 

For this model, n = 0 and the average scale factor a = /c 2 exp(//) turns the 
metric coefficients A, B and C into 


A 

B 

C 


Pik 2 exp(//) exp 
^ 2^2 exp(//) exp 
P 3 /C 2 exp(//) exp 


qi exp{—3lt) 

3/^7 J’ 

q' 2 exp(- 3 //) 

3/fc? J’ 

g 3 exp(- 3 /t)' 

3lk2^ 


(50) 

(51) 

(52) 
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Figure 2: Behavior of w versus time for f > 0 with n = 2, A = 1, Z = 1, /ci = 
0, qi = l = q 2 and ga = -2. 


The directional Hubble parameters Hi become 


Qi 


Hi = l + -^exp(—3/t). 
C2 


(63) 


The mean generalized Hubble parameter and volume scale factor turn out to 
be 

H = l V = k2^exp{3lt). (54) 

The mean anisotropy parameter, expansion scalar and shear scalar are 


A = 


qi^ + q2^ + 


e = 31, = 


gi" + q2^ + gs" 


exp(6Zt) ’ ’ 2 /c 2 ^ exp(6/t) 

The energy density and pressure of the universe take the form 


1 


P = 


12(A + 27r)(A + 47r) _ 

gi^ + g2^ + ga 


4(A + 37r)<^3Z^ + 


qiq2 + q2q3 + qsqi 

k 2 ^ exp(6Zf) 


- A<^3Z^ + 


/c 2 ° exp(6/f) 


P = 


-i 

12(A + 27r)(A + 47r) 




4/ex 


^ fc2®exp(6/f) 


:p(6Zf) 


(55) 


(56) 


(57) 
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Figure 3: Behavior of energy density and pressure versus time for t > 0 with 
n = 2, A = 1, / = 1, ^2 = 1, qi = I = q 2 and qs = -2. 



Figure 4: Behavior of w versus time for t > 0 with n = 2, A = 1, I = 
1, k 2 = l, qi = l = q 2 and qs = -2. 

For this model, the plots of p, P and w against time coordinate t are shown 
in hgure 3 and 4 respectively. It can be seen from hgure 4 that w —>■ — 1 as 
f —>■ cx) which indicates that the non-singular model corresponds to a vacuum 
fluid dominated universe. 


4 Concluding Remarks 

This paper is devoted to discuss the current phenomenon of accelerated ex¬ 
pansion of universe in the framework of newly proposed f{R,T) theory of 
gravity. For this purpose, we take f{R, T) = R + 2AT and explore the exact 
solutions of Bianchi type I cosmological model. We obtain two exact solu- 
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tions using the assumption of constant value of deceleration parameter and 
the law of variation of Hubble parameter. The obtained solutions correspond 
to two different models of universe. The hrst solution forms a singular model 
with power law expansion while the second solution gives a non-singular 
model with exponential expansion of universe. The physical parameters for 
both of these models are discussed below. 

The singular model of the universe corresponds to u 7 ^ 0 with aver¬ 
age scale factor a = {nit + ki)«. This model possesses a point singularity 
when t = tg = The volume scale factor and the metric coefficients 

A, B and C vanish at this singularity point. The cosmological parameters 
Hi, H2, H3, H, 6, and are all inhnite at this point of singularity. If we 
choose ki = 0 , hgure 1 suggests that energy density of the universe is zero at 
this time. The pressure approaches negative inhnity as t —)■ 0. This strong 
negative pressure is an indication of dark energy. For this model, tc —| 
as t —)■ 00 which corresponds to a radiation dominated universe. The mean 

anisotropy parameter A also becomes inhnite at this point for 0 < n < 3 and 

2 

vanishes for n > 3. Moreover, the isotropy condition, i.e., ^—?• 0 as f —)■ cx), 
is verihed for this model. All these conclusive observations suggest that the 
universe starts its expansion with zero volume, strong negative pressure and 
energy density from t = tg and it will continue to expand for 0 < n < 3. 

Now we discuss the non-singular model of the universe corresponds to 
u = 0. For this model the average scale factor is a = A: 2 exp(Zt). The non¬ 
singularity is due to the exponential behavior of the model. The expansion 
scalar 9 and mean generalized Hubble parameter H are constant in this 
case. For hnite values of t, the physical parameters Hi, H2, H^, cx^ and 
A are all hnite. The metric functions are dehned for hnite time and the 
isotropy condition is satished. There is an exponential increase in the volume 
as the time grows. However, energy density is approximately zero initially 
and becomes constant after some time. Pressure of the universe remains 
in the negative zone for this model which may be an indication of presence 
of dark energy in the universe. Figure 4 suggests that tc —)■ — 1 as f —>■ 
00 . Thus the exponential model corresponds to a vacuum huid dominated 
universe. According to the observations 0 . the expansion of the universe 
is accelerating when w — 1 . 

Therefore, it is hoped that the problematic issues such as dark energy and 
accelerated expansion of universe may be addressed using modihed theories 
of gravity especially f{R, T) gravity. It would be interesting to explore more 
Bianchi type solutions in this context. Exact solutions of Bianchi type V 
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cosmological model in this theory are under process. 


Acknowledgement 

The author is thankful to National University of Computer and Emerging 
Sciences (NUCES) Lahore Campus for funding the PhD programme. The 
author is also grateful to the anonymous reviewer for valuable comments and 
suggestions to improve the paper. 


References 

[1] Harko, T., Lobo, F.S.N., Nojiri, S. and Odintsov, S.D.; Phys. Rev. 
D84(2011)024020. 

[2] Riess, A.G. et ah: Astron. J. 116(1998)1009; Riess, A.G. et ah: Astro- 
phys. J. 607(2004)665. 

[3] Perlmutter, S. et ah: Astrophys. J. 517(1999)565. 

[4] Carmeli, M.: Commun. Theor. Phys. 5(1996)159. 

[5] Seprgel, D. N. et ah: Astrophys. J. Suppl. 148(2003)175; Spergel, D.N. 
et ah: Astrophys. J. Suppl. 170(2007)377. 

[6] Allen, S. W. et ah: Mon. Not. R. Astron. Soc. 353(2004)457. 

[7] Bennett, C. L. et ah: Astrophys. J. 148(2003)1. 

[8] Hogan, J.: Nature 448(2007)240. 

[9] Tegmark, M. et ah: Phys. Rev. 069(2004)103501 

[10] Abazajian, K. et ah: Astron. J. 129(2005)1755. 

[11] Astier, P. et ah: Astron. Astrophys. 447(2006)31. 

[12] Nojiri, S. and Odintsov, S.D.: Int. J. Geom. Meth. Mod. Phys. 
4(2007)115. 


14 



[13] Turner, M.S., Huterer, D.: J. Phys. Soc. Jap. 76(2007)111015; Frie- 
man, J., Turner, M. and Huterer, D.: Ann. Rev. Astron. Astrophys. 
46(2008)385. 

[14] Li, M., Li, X.D., Wang, S. and Wang, Y.: Commun. Theor. Phys. 
56(2011)525. 

[15] Copeland, E.J., Sami, M. and Tsujikawa, S.: Int. J. Mod. Phys. 
D15(2006)1753. 

[16] Sahni, V. and Starobinsky, A.: Int. J. Mod. Phys. D9 (2000) 373; Sahni, 
V.: Lect. Notes. Phys. 653 (2004) 141. 

[17] Carroll, S.M.: Living Rev. Rel. 4(2001)1. 

[18] Weinberg, D.H.: New. Astron. Rev. 49(2005)337. 

[19] Straumann, N.: Mod. Phys. Lett. A21 (2006) 1083. 

[20] Yang, R.J.: Europhys. Lett. 93(2011)60001; Wei, H., Ma, X.P. and 
Qi, H.Y.: Phys. Lett. B703(2011)74; Wu, P.X. and Yu, H.W.: 
Eur. Phys. J. 071(2011)1552; Wu, P.X. and Yu, H.W.: Phys. Lett. 
B703(2011)223; Bamba, K., Ceng, C.Q., Lee, C.C and Luo, L.W.: 
JCAP 1101(2011)021; Li, B., Sotiriou, T.P. and Barrow, J.D.; Phys. 
Rev. 083(2011)104017. 

[21] Felice, A.D and Tsujikawa, S.: Living Rev. Rel. 13(2010)3; Sotiriou, 
T.P. and Faraoni, V.: Rev. Mod. Phys. 82(2010)451; Clifton, T., Fer¬ 
reira, P.G., Padilla, A. and Skordis, C.: Phys. Rept. 513 (2012)1; 
Nojiri, S. and Odintsov, S.D.: Phys. Rept. 505(2011)59. Bamba, K., 
Capozziello, S., Nojiri, S. and Odintsov, S.D.: Astrophys. Space Sci. 
342(2012)155. 

[22] Sharif, M. and Zubair, M.: Adv. High Energy Phys. 2013(2013)790967. 

[23] Sharif, M. and Kausar, H.R.: JCAP 07(2011)022. 

[24] Sharif, M. and Kausar, H.R.: Int. J. Mod. Phys. D20(2011)2239. 

[25] Sharif, M. and Kausar, H.R.: Astrophys. Space Sci. 331(2011)281. 

[26] Sharif, M. and Kausar, H.R.: Mod. Phys. Lett. A25(2010)3299. 


15 



[27] Bamba, K., Nojiri, S., Odintsov, S.D. and Saez-Gomez, D.: Phys. Lett. 
B730(2014)136. 

[28] Bamba, K., Makarenko, A.N., Myagky, A.N., Nojiri, S. and Odintsov, 

S. D.: JCAP 01(2014)008. 

[29] Bamba, K., Nojiri, S. and Odintsov, S.D.: Phys. Lett. B698(2011)451. 

[30] Capozziello, S. and Vignolo, S.: Int. J. Geom. Meth. Mod. Phys. 
8(2011)167. 

[31] Gapozziello, S., Darabi, F. and Vernieri, D.: Mod. Phys. Lett. 
A26(2011)65. 

[32] Gapozziello, S., Laurentis, M.D., Odintsov, S.D. and Stabile, A.: Phys. 
Rev. 083(2011)064004. 

[33] Elizalde, E., Nojiri, S., Odintsov, S.D. and Saez-Gomez, D.: Eur. Phys. 
J. 070(2010)351. 

[34] Bamba, K., Geng, G., Nojiri, S. and Odintsov, S.D.: Mod. Phys. Lett. 
A25(2010)900. 

[35] Gapozziello, S., Laurentis, M.D., Nojiri, S. and Odintsov, S.D.: Gen. 
Rel. Grav. 41(2009)2313. 

[36] Multamaki, T. and Vilja, L: Phys. Rev. 074(2006)064022; Multamciki, 

T. and Vilja, L; Phys. Rev. 076(2007)064021. 

[37] Gapozziello, S., Stabile, A. and Troisi, A.: Glass. Quantum Grav. 
24(2007)2153. 

[38] Hollenstein, L. and Lobo, F.S.N.: Phys. Rev. 078(2008)124007; Sharif, 
M. and Kausar, H.R.: J. Phys. Soc. Jpn. 80(2011)044004; Shojai, A. 
and Shojai, F.: Gen. Relativ. Gravit. 44(2011)211. 

[39] Azadi, A., Momeni, D. and Nouri-Zonoz, M.: Phys. Lett. 

B670(2008)210; Momeni, D. and Gholizade, H.: Int. J. Mod. Phys. 
018(2009)1719; Sharif, M. and Arif, S.: Astrophys. Space Sci. 
342(2012)237. 

[40] Sharif, M. and Shamir, M.F.: Mod. Phys. Lett. A25(2010)1281. 


16 



[41] Sharif, M. and Shamir, M.F.: Class. Quantum Grav. 26(2009)235020; 
Sharif, M. and Shamir, M.F.: Gen. Relativ. Gravit. 42(2010)2643. 

[42] Shamir, M.F., Jhangeer, A. and Bhatti, A.A.: Chin. Phys. Lett. 
29(8) (2012)080402. 

[43] Myrzakulov, R.: Phys. Rev. D84(2011)024020. 

[44] Sharif, M. and Zubair, M.: JCAP 03(2012)028. 

[45] Sharif, M. and Zubair, M.: J. Phys. Soc. Jpn. 82(2013)064001. 

[46] Houndjo, M.J.S.: Int. J. Mod. Phys. 021(2012)1250003. 

[47] Adhav, K.S.; Astrophys. Space Sci. 339(2012)365. 

[48] Reddy, D.R.K., Santikumar, R. and Naidu, R.L.: Astrophys. Space Sci. 
342(2012)249. 

[49] Reddy, D.R.K., Santikumar, R. and Pradeepkumar, T.V.: Int. J. Theor. 
Phys. 52(2013)239. 

[50] Naidu, R.L., Reddy, D.R.K., Ramprasad, T. and Ramana, K.V.: Astro¬ 
phys. Space Sci. 348(2013)247. 

[51] Ahmed, N. and Pradhan, A.: Int. J. Theor. Phys. 53(2014)289. 

[52] Landau, L.D. and Lifshitz E.M.: The Classical Theory of Fileds 
(Butterworth-Heinemann, 2002). 

[53] Berman, M.S.: Nuov. Cim. B74(1983)182; Shamir, M.F. and Bhatti, 
A.A.: Can. J. Phys. 90(2012)193; Singh, C.P. and Kumar, S.: Int. J. 
Theor. Phys. 47(2008)3171-3179; Singh, C.P., Zeyauddin, M. and Ram, 
S.: Int. J. Theor. Phys. 47(2008)3162-3170. 

[54] Hogan, J.: Nature 448(2007)240. 


17 



